Suppression of Tunneling of Superconducting Vortices Caused by a Remote Gate: 

Example of an Extended Object Tunneling 
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We discuss a recent experiment in which the resistance of a superconducting film has been mea- 
sured in magnetic field. A strong decrease of the superconducting film resistance has been observed 
when a metallic gate is placed above the film. We study how the magnetic coupling between vortices 
in a thin superconducting film and electrons in a remote unbiased gate suppresses the tunneling rate 
of the vortices. We examine two general approaches to analyze tunneling in the presence of slow 
low-energy degrees of freedom: the functional-integral and scattering formalisms. In the first one, 
the response of the electrons inside the metallic gate to a change in the vortex position is described 
by the "tunneling with dissipation". We consider the Eddy current induced in the gate by the 
magnetic flux of the vortex as a result of tunneling. In the second approach, the response is given 
in terms of scattering of the electrons by the magnetic flux of the vortex in a way similar to the 
Aharonov-Bohm scattering. A sudden change in the vortex position leads to the Orthogonality 
Catastrophe that opposes the vortex tunneling. We show that the magnetic coupling between the 
vortices and the electrons inside the gate can lead to a dramatic suppression of the vortex tunneling, 
restoring the superconducting property in accord with the experiment. 



I. INTRODUCTION 



The vortex motion in superconductors is a source 
of energy losses, destroying the perfect conductivity 
of the superconductors. The dissipation is caused by 
non-superconducting electrons located inside the vortex 
core. 1 The pinning potential created by impurities op- 
poses the motion of vortices. This potential results from 
the action of the impurities on the vortex core averaged 
over the area of the core. Therefore, the pinning potential 
has minima typically separated by a distance of order of 
the coherence length £. 2 ' 3 Vortices may change their po- 
sitions either by thermal activation 4 or by quantum tun- 
neling between the potential minima at low enough tem- 
peratures. 5-8 In this context, the observation of a strong 
decrease of the resistance when an unbiased metallic gate 
is placed above an amorphous superconducting film 9 is of 
great interest. We believe that this experiment provides 
a strong argument that at low temperature the motion of 
vortices is indeed realized by quantum tunneling (in the 
experiment described in Ref. 9 it occurs at T < O.IK). 
If so, the tunneling of a vortex in a thin superconducting 
film is a unique example of tunneling of a very extended 
"object". 

In the discussed experiment, the resistance of a su- 
perconducting film has been measured at various mag- 
netic fields, both with and without a gate. In the ab- 
sence of a gate, in magnetic fields lower than the crit- 
ical one (H < Hci) > the resistance initially decreases 
with lowering the temperature, but eventually saturates 
at a finite value. The saturation indicates the possibil- 
ity of vortex tunneling. When an unbiased metallic gate 
is placed above the superconducting film the resistance 
reduces significantly with no indication of saturating at 
a finite resistance when T — > 0. Remarkably, the effect 
of the gate becomes noticeable at the same temperatures 
where the resistance of the ungated film starts to satu- 



rate. It is worth mentioning that the gate is separated 
from the film by an oxide layer of 160A Therefore, the 
film is thermally isolated from the gate ruling out the 
possibility that the saturation of the resistance in the 
ungated film can be attributed to heating. In this pa- 
per we identify the mechanism causing the suppression 
of the vortex motion in the presence of the gate, which is 
effective in the tunneling regime only. The fact that the 
finite resistance at low temperatures has been eliminated 
by placing a remote isolated gate strongly confirms the 
tunneling nature of the vortex motion. 

In the experiment of Ref. 9 the film thickness is a w 
30A and the gate thickness is d s=s 400A An important 
feature indicating that the gate and the film are well 
separated, is that the superconducting transition tem- 
perature, T c , as well as the critical magnetic field, are 
practically unchanged by adding the gate. Since the gate 
does not affect the superconductivity at T ~ T c when 
it is the weakest, its influence on the superconducting 
properties, like the energy gap, at lower temperatures 
can be ignored. One should also keep in mind that the 
gate does not influence the resistance of the film when the 
vortex motion is thermally activated. In view of all of the 
above, we concentrate only on examining the influence of 
the gate on the vortex tunneling rate. We assumed that 
the superconducting film and the gate are magnetically 
coupled via the magnetic field of the vortices that pierces 
through the gate. 

We employ here the following strategy. We accept the 
tunneling of vortices at low temperatures as an estab- 
lished experimental fact. We do not try to calculate the 
tunneling rate of the vortices. Instead, we concentrate 
on how the response of the electrons inside the gate to 
a change of the vortex position suppresses the tunneling 
rate. With this question in mind, in the complex prob- 
lem of the vortex tunneling, we wish to isolate the effect 
induced by the gate. 

In fact, little is known about the motion of vortices at 
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FIG. 1: The vortex in the effective potential landscape is rep- 
resented by a hopping "particle", (a) The tunneling rate in the 
absence of the gate Y un G exceeds the energy difference between 
the potential minima s; the vortices are in a "metallic" phase 
(b) The tunneling rate is reduced by the gate to Yg < s, and 
the system of vortices becomes an "insulator". 



low temperatures. Fortunately, for studying the role of 
the gate, it is sufficient to assume that the change in the 
vortex position is a discrete tunneling event. This can 
be a tunneling of a single vortex, a bundle of vortices, or 
topological defects such as dislocation pairs in the case 
of a vortex lattice (or a glass state) . Phcnomenologically, 
the change in the vortex position can be described by the 
hopping Hamiltonian 

H = Siafai + U '<..",' "./ + h - c ) ■ (!) 

i (hi) 

According to the standard criterion of the Metal- 
Insulator transition 11 , the ratio of the variation of the 
potential minima e = (s^ to the typical value of the tun- 
neling rates T = (Tij) specifies whether the vortices are 
itinerant or localized. The finite resistance at low tem- 
perature in the absence of the gate indicates that the 
vortices are mobile, i.e., the tunneling rate in the un- 
gated film T U71 q > e. The dissipationlcss nature of the 
superconducting film is revived when the tunneling rate 
is reduced by the gate to Tq < £', see Fig. 1 which illus- 
trates the two cases. Thus, we interpret the experiment 9 
as a transition from "metallic" to "insulating" phases in 
a system of tunneling vortices that is induced by the gate. 
(Because of the strong interaction between the vortices, 
it may be enough to have a fraction of vortices with a 
suppressed probability of tunneling acting as the pinning 
centers for the entire ensemble of vortices.) 

The fact that the tunneling of the vortices can be 
blocked by placing a gate above the film indicates that 
the tunneling event gives rise to a dramatic response of 
the electrons inside the gate. This response can be an- 
alyzed in terms of low-energy electron-hole excitations 
"decorating" the tunneling event. One may consider the 
cloud of virtual excitations as part of the tunneling pro- 
cess that lasts long after the change in the vortex posi- 
tions occurs. In the present work we concentrate on the 
response of the electrons to the tunneling that can lead 
to the strong suppression of the tunneling rate. 

One may better understand the specifics of tunneling 
with the participation of low-energy degrees of freedom, 
i.e., in the presence of dissipative environment, using the 
picture given by Iordanskii and one of the authors. 12 
Originally it described the quantum formation of a nu- 



clcation center in the decay of a metastable macroscopic 
state, but it could also provide a general perspective. 
The quantum nucleation is an example of a tunneling 
process in which a large number of degrees of freedom 
participate. Alternatively, this kind of process can be 
treated as a tunneling of an artificial "particle" in a multi- 
dimensional space. When low-energy degrees of freedom 
are involved in the process of quantum nucleation, the 
nucleation develops in two stages. 12 Namely, the motion 
of the " particle" in this multi-dimensional space along the 
trajectory minimizing the imaginary time action consists 
of fast and slow stages. 

The slow stage appears because of the long time needed 
for the slow low-energy degrees of freedom to adjust 
themselves to the new state of the fast degrees of free- 
dom. 13 This time is much longer than needed for fast 
the degrees of freedom to complete the tunneling. That 
is why the tunneling has to develop in two stages. It 
has been shown in Rcf. 12 that despite the fact that the 
low-energy degrees of freedom yield only a small contri- 
bution to the energy of the barrier, their participation in 
the tunneling process increases parametrically the over- 
all tunneling time. This results in a large increase of the 
action and, correspondingly, in the strong suppression of 
the tunneling rate. 14 

The described picture of changing a quantum state in 
the presence of low-energy degrees of freedom is rather 
typical for condensed matter systems. In the course of 
the fast stage of the process a quantum mechanical ob- 
ject changes its state (position, spin projection, phase of 
the Josephson junction, etc). The accompanying slow 
degrees of freedom act as an environment for the fast 
degrees of freedom. In the discussed problem of tunnel- 
ing in the gated superconducting film, the tunneling of 
the vortex from one potential minimum to another cor- 
responds to the fast stage. During the slow stage, the 
electrons inside the gate adjust their state to the new 
position of the vortex. 

The ensemble of electron-hole pairs in the gate repre- 
sents the low-energy degrees of freedom of the environ- 
ment. The environment produces the most significant 
effect on tunneling at large time differences when the 
slow degrees of freedom have enough time to develop. 
Naturally, the effect is the strongest when the tunnel- 
ing occurs between states that are almost degenerate. In 
the latter case the long time response can considerably 
reduce or even block the tunneling 15,16 (this statement 
is often formulated in terms of the dissipative Quantum 
Phase Transition 17 ). With this in mind, we concentrate 
only on the slow stage that develops when the tunneling 
of the vortex degrees of freedom is mostly accomplished 
without specifying how the fast stage develops. 

In this paper, we examine two different approaches to 
describe the response of the electrons inside the gate to 
the change of the vortex position. In Sees. II and III 
we consider the dissipative Eddy (Foucault) currents in- 
duced in the gate which continue long after the vortex 
changed its position. We formulated the effect of the 



gate on the vortex tunneling in terms of the effective 
action of a vortex. We integrate out the environmen- 
tal degrees of freedom inside the gate, and obtain the 
dissipative term in the effective action. This allows us 
to consider the vortex tunneling in the gated supercon- 
ducting film in the context of the well-known problem of 
tunneling in the presence of a dissipative environment 18 , 
Sec. IV. Alternatively, one can analyze the slow stage 
in terms of the scattering of electrons. In Sec. V we de- 
scribe the elastic scattering of the gate electrons on the 
vector potential of the magnetic field of the vortex. The 
zero overlap between the states of the electrons before 
and after the change of the scattering potential is known 
generally as the Orthogonality Catastrophe 19 (OC). We 
show that the OC caused by the change in the vortex 
position, effectively suppresses the vortex tunneling. The 
novel element here 20 is that the OC is connected to the 
Aharonov-Bohm effect. 21,22 We obtained that the OC is 
significantly more effective in suppressing the tunneling 
rate than the Eddy current. In Sec. VI we discuss the re- 
lation between the two approaches and the peculiarities 
of the tunneling of an extended object, such as a vortex 
in a thin superconducting film. 



II. EDDY CURRENTS INSIDE THE GATE 

The magnetic field of the vortex inside a supercon- 
ducting film is similar to the magnetic field of a solenoid 
with a radius equal to the magnetic penetration depth A 
which can be very large. 23 Outside of the film, the mag- 
netic field decays as a function of the height and deflects 
into the radial direction 24 ' 25 : 



A. vo) -(r,z, R.) 



iqr-q\z\ fol X g 



(27T) 2 



q 2 (\- 1 + 2 q y 



(2) 

Here r is the radial vector in cylindrical coordinates with 
the origin at the vortex center R. Since a vortex in a 
film can move only in the x — y plane, R is a two di- 
mensional vector. The parameter a = 1/2 is the total 
flux of a vortex in a superconductor measured in units 
$ = 2-irhc/e. 

In the case of a thin superconducting film, a Pearl vor- 
tex is a macroscopically large object. The penetration 
depth A is 23 : 



A = A 2 £,/a A ;» a, 



(3) 



where a is the thicknesses of the superconducting film, 
X-3D is the penetration depth in a disordered bulk su- 
perconductor. In Eq. (3), A3D ~ ^Liis/l) 1 / 2 , where 
Al = (mc 2 /47rne 2 ) 1 / 2 and £0 are the London length and 
the coherence length of a clean superconductor; I is the 
mean free path. 

In the experiment of Ref. 9 the gate thickness d is 
much smaller than A. Because of the exponential decay 
of the Fourier components in the z-direction, see Eq. (2), 
momenta that contribute mostly are limited to q < 1/d. 




FIG. 2: Superconducting film magnetically coupled to a metal- 
lic gate. The magnetic field of the vortex pierces through the 
gate. In the experiment of Ref. 9 d ~ 400A, a » 30A, and the 
insulating layer between the gate and the film is 160A thick. 



More accurately, q < min{d _1 , Sr -1 }, where Sr is the 
typical distance that a vortex has to tunnel, Sr ~ £. 
To avoid unnecessary complications, we ignore the space 
between the superconducting film and the gate since it is 
considerably smaller than the gate thickness, see Fig. 2. 
Then, for small momenta that we are interested in, the 
deflection of the magnetic field in the space between the 
film and the gate can be neglected. 

To describe the tunneling of a vortex in the presence 
of a gate one has to deal with an imaginary time action: 



S — S sc + S„ate + $11 



(4) 



In what follows we discuss each term in the action S. 

The term S sc is an action of the superconducting film 
in the absence of the gate. Since we are not trying to 
solve the problem of the vortex tunneling in full scale, 
but are interested only in the effect of the gate on the 
tunneling rate, this part of the action is not specified. 

In writing S ga te which describes the dynamics of the 
electrons in the gate one should keep in mind the follow- 
ing argument. The charge and current densities relevant 
for the long tail response of electrons inside the gate to 
the tunneling of the vortex are characterized by large 
length and time scales. Therefore, their dynamics can 
be described macroscopically. Since the deviations of the 
charge and current densities from their equilibrium values 
are small, the action that describes the dynamics of their 
fluctuations should be consistent with the Fluctuation- 
Dissipation Theorem (FDT). (Examples of such an ap- 
proach can be found in the calculation of the dephasing 
time of the cooperons due to the electromagnetic fluctua- 
tions 26 , and also in a macroscopic calculation of the zero 
bias anomaly. 27 ) 

The current in the gate has two contributions. One is 
the Ohmic response to the electric field, while the other 
one is the diffusive current from the gradients of the den- 
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sity; J(r,r) = 3 hmic(r,r) - DVp(r,r). The fluctua- 
tions of the charge and current densities can be expressed 
through the correlation function of the Ohmic part of the 
current written in terms of the Matsubara frequency as 
follows: 

Kij (k, iu n ) = (Jofcmic( k > iuj n)3i hmlc (-k, -iu n )) (5) 
= crjj(k, oj n )\oj n \ + <7i,i'(k, LJ n )D jtj >(\s.,oj n )ki>kj>. 

Here the diffusion constants tensor D and the conduc- 
tivity a are connected through the Einstein's relation 
(7 = e 2 (dn/dp)D. The gate, being a simple homoge- 
neous metal, is adequately described by the Drude for- 
mula. We assume that the external magnetic field is clas- 
sically weak, u> c t C 1, so that we can ignore the Hall 
conductivity in the gate; r is the mean free time in the 
gate. Since we are interested in the low frequency and 
long wave length behavior only, we take the conductivity 
a to be constant. 

Following the above arguments, one can write the gate 
action S ga t e for the charge and current densities inside 
the gate in the Matsubara time: 

Sgate^^ J dTidr 2 j dr 1 dr 2 L FDT [r 1 ,Ti;r 2 ,T 2 } (6) 

+ J dr J dr[C 

cont ax J 

The first term in the above action describes the charge 
and current density fluctuations in accordance with the 
FDT: 

Lfdt = J(ri,n) +DVp(r 1 ,Ti) Kfa -r 2) r 1 -r 2 ) 

(7) 

J(r 2 ,T 2 )+Wp(r 3 ,r 2 ) . 

Since the fluctuations of the current and density are not 
independent, the second term in Eq. (6) imposes the 
charge continuity condition with the use of the Lagrange 
multiplier (f>: 



r,r 



(ip(r,r)+VJ(r, r)) 



(8) 



The term CMaxweii describes the interactions of the elec- 
tromagnetic fields and the charge and current densities in 
the gate in a way that reproduces the Maxwell equations: 



C-Max — i— [JA-ind + Pfind] 
C 



(9) 



+ 



8?r 



(V x A md f 



The factor i in the coupling of the electromagnetic fields 
with the charge and current densities appears because of 
the imaginary time. 28 The plus in the magnetic field term 
(which is a consequence of the i-factor) is needed to get 
the repulsive sign in the Ampere interaction of currents, 



like in the Coulomb interaction of charges. As usually, 
the total potentials are the sum of the external and in- 
duced potentials; in the discussed system the external 
source is the field of the vortices. 

The last term in the action (4), Si n t, describes the con- 
nection between the superconducting film and the gate. 
The current and charge densities in the gate interact with 
the vector and scalar potentials created by the supercon- 
ducting film: 

Sm = i J dr y dr |ij(r, r) • A sc (r, r) 

+ p(r,r) v?sc (r,r)|(10) 

We are interested in the limited problem of the long tail 
response of electrons inside the gate that develops when 
the tunneling of the vortex degrees of freedom is mostly 
accomplished. Therefore, we consider large enough r 
when one can assume that the deformation of the field 
A sc (r, r) appearing during the process of vortex tunnel- 
ing has been already relaxed. [By deformation we mean 
the deviation of A sc (r, r) from the field of the "rigid" 
vortex centered at R(r), as it is given by Eq. (2). Here 
R(t) denotes the position of the vortex at time r.] With 
this in mind, we put in Si n t 



A sc (r,r) = A„ or (r;R(r)). 



(11) 



Furthermore, we ignore the scalar potential, y> sc (r, t) = 
ip vor (r; R(t)) = 0, relying on the known fact that the re- 
distribution of the charge density around the vortex core 
is negligible. We still have to justify our treatment of the 
electromagnetic field A sc as a given external field in the 
analysis of the long tail response of the gate electrons. 
We will see in the next section that the magnetic field 
created by the low-frequency components of the dissipa- 
tive Eddy currents is much smaller than the field of the 
vortex, and can perturb the superconducting film only 
weakly. Therefore, the gate does not provide a substan- 
tial feedback effect to the superconducting film during 
the slow stage. 

To get the response of the environment on the vortex 
motion, one has to integrate out the gate degrees of free- 
dom. Since S ga t e and Si n t are quadratic in the charge 
and current densities, this immediately results in 

Senv = ^ J dr J dr^-3 c i(r,T) ■ A„ or (r, R(t)). (12) 

The current J c ;(r, r) has to be found by solving the clas- 
sical equations of motion. Since the current J(r, r) de- 
scribes the long-time response of the electrons in the gate, 
the term S env is non-local in time. 



III. THE SOLUTION OF THE EQUATION OF 
MOTION 



As we stated above, the main contribution to the in- 
duced current emerges from the components of the field 
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with q < 1/d for which the depletion of the magnetic field 
of the vortex from the z direction is negligible. Therefore, 
we can ignore the current component in the z-direction 
and consider only the components that are parallel to 
the plane. With this in mind, we decompose the current 
parallel to the plane of the gate into two components: a 
longitudinal component along the two dimensional vec- 
tor q and a transverse component perpendicular to it, 
J(q, t) = j"(q, r)q + J ± (q, t)z x q. In the following 
we use the Fourier transformation for the in-plane coor- 
dinates only and keep the vertical coordinate separately. 
Then, the kernel in Eq. (5) is diagonal in the chosen 
basis. As it follows from Eq. (2), the vector potential 
A vor (r) contains the transverse component only. There- 
fore, since the use of the FDT for describing the slow 
stage implies a linear response to the external field, the 
longitudinal component j" cannot be generated. This 
is a direct consequence of the magnetic coupling between 
the superconducting film and the gate. It is very different 
from the zero bias anomaly in which only a longitudinal 
current is generated. 27 

The variation of S ga t e + Si n t with respect to the current 
and charge densities determines the equations of motion. 
We use these equations for finding the transverse current: 

J dz'K^iq; z, z'; iw n ) J x (q, z', iuj n ) (13) 

If 1 

= j / dz'U{q, z-z')J J '(q, z',iu n ) - i-A^ or (q,iuj n ). 

c z J c 

The L.H.S of this equation is the total vector potential 
in the transverse direction calculated via its Ohmic re- 
sponse: J = (auo n / c)A^ otal . The first term in the R.H.S 
is the vector potential of the induced field and the second 
one is the external potential of the vortex. Neglecting rel- 
ativistic effects, the instant kernel U can be written as 
the Fourier transform of the instant in time Biot-Savart 
kernel U(r — r') = l/|r — r'| with respect to the in-plane 
momenta, £/(q, z — z') = 2ne~ q \ z ~ z ' jq. 

The above equation for the transverse current can be 
rewritten in the form: 

f I I 

/ dz'i ± (q;z,z';iw„)J ± (q, iu n ) = -i — A x r (q, iu n ), 

(14) 

where the kernel L is 

L ± (q;z,z';i^n) = (15) 
f 5{z - z') + 2 -H^A e - q \z-z'\/ q o < z, z' < d; 
\ otherwise. 

To find the transverse current one has to invert this ker- 
nel. 

In analogy to the skin-effect, one can define a screen- 
ing length <5(q, u n ) = 1/^q 2 + Ai:a\Lo n \/ c 2 . [Usually, the 
skin-effect is discussed in the case of an electromagnetic 
wave propagating normally to a surface of a metallic slab. 
The geometry of the problem studied here is different as 
the magnetic field is normal to the slab surface, while 



the propagation is parallel to it. Still, the surface current 
appearing in the gate screens out the high-frequency com- 
ponents of the field in the bulk of the slab.] The kernel 
L 1 - can be inverted in the two limits: (i) for the com- 
ponents of the electromagnetic field that are transparent 
for the gate, d/<5(q, io n ) <C 1; or (ii) for the components 
with <i/<5(q, oj n ) 3> 1 that are well screened by the sur- 
face currents. These are the thin and thick gate limits, 
respectively. 

In the thin gate limit, the change in the current along 
the z direction is minor, and a reasonable approximation 
(up to liner terms in qd) is to consider the current to 
be homogenous in the z direction, I-id = dJ^(q, z, iui„). 
Then, 

l2D = ~* 2™™\u I — A vor(^ iuJ n) + 0{qd). (16) 

1+ ° 

One can observe that the expression for the current I2D 
is identical to the current in a two-dimensional system 
with a 2D = ad, and where 2na 2D \uj n \/qc 2 is the cur- 
rent screening operator. [The term "current screening 
operator" is used here in analogy with the polarization 
operator. It describes the screening of the transverse 
component of the time dependent vector potential by the 
induced currents.] 

At low frequencies such that 2Tra 2D \uj n \/c 2 < q < rf" 1 
(this automatically implies the thin gate limit), the cur- 
rent in the gate screens weakly the field produced by the 
vortex. The total field is approximately just the field of 
the vortex, and the current is merely the Ohmic response 
to it: 

I 2D = -i^^ or (q, w n ). (17) 
c 

For higher frequencies, 2ira 2D \uj n \/ c 2 > q, but still in the 
thin gate limit, the situation is rather different. Since 
the effect of dissipation is stronger for better conducting 
gates, we are interested in the case when 2-ku 2D jc S> 1. 
Then, there is a window qc 2 /2-ko~ 2D < \uj n \ < qc in which 
the current I2D obeys a London-like equation 24 , I2D = 
~ *ff Av 0r (<l, i^n). Since the electrons instantly respond 
to the potential, the contribution to the action in this 
limit is not of a dissipative character. Rather, it provides 
a local in time term acting as an additional potential that 
should be added to S sc . This kind of contribution is not 
considered here. 

For a thick gate, i5(q, ui n ) <IC d, the limit of low fre- 
quencies does not exist for q < d . In this limit the 
screening length is equal to 5(tu n ) = cj \J 'A-na\uj n \. At 
such high frequencies the current flows in the reduced 
volume (which effectively is a thin slab of width 8(aj n )) 
as follows: 

J x (q,^ia;„) = -i^ 1 ^ Lo n )e~ z l s ^)+^ 

x^ i Dr (q,z,^„) + 0(e- d/5 ) ■ (18) 
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Notice that the Fourier components A^ or (q, z) decay ex- 
ponentially on z as e~ qz . Therefore, the factor e qz in 
the solution above is canceled out leaving the induced 
surface current with the z-dependence e~ z ^^ q, ^ n \ Inte- 
grating Eq. (18) in the z-direction yields the London- like 
current I2D = —i^-^vor( c l^ u, n) exactly as in the thin 
gate limit. 

So far, we have ignored relativistic effects that ap- 
pear at high frequencies when \uj n \/c > q. In this 
case, one has to substitute q by the relativistic com- 
bination \J q 2 + (uj n /c) 2 . Using the relativistic equa- 
tions of motion, one can show that the current J c ; is 
like in the Ohmic regime but with a 2D replaced by 
a 2D /{\ + 2na 2D /c). In the limit 2na 2D /c > 1, the 
dissipation caused by the Cherenkov's radiation 29 corre- 
sponds to the effective conductivity equal to c/2-k. Still, 



the effect of the relativistic region is negligible in com- 
parison to the low- frequency contribution to S env . 

To conclude, let us come back to our assumption about 
the absence of a feedback from the gate to the super- 
conducting film. As we have showed above, the most 
significant contribution to S env originates from the re- 
gion 2ir\u> n \(j 2D I (gc 2 )<Sl when the Eddy current is in 
the Ohmic regime. In this limit, the low-frequency com- 
ponents of the vortex field are poorly screened by the 
induced current in the gate. Hence, the feedback from 
the gate to the superconducting film can be neglected. 

We are ready to obtain the dissipative term in the ac- 
tion describing the effective response of the environmen- 
tal degrees of freedom on the vortex tunneling. Inserting 
the current J x from Eq. (17) into Eq. (12), one gets 



Senv = / ' dTldT2 J2\ 



-iU] n (T 1 -T 2 ) 



d2q .rf ze tq(R(Ti) - R(T2)) ~ iq x A ™ r ("~ q ' z ) iq x A ™^(q, 



(27T? 



(19) 



Performing the sum over the frequencies and rewriting the action in terms of the magnetic field we get eventually 

d 2 q J e ii(R(^)-R(r 2 )) Bi or (-<i,z)Bi or (q,z) 



7T<7 

' 2/3 2 c 2 



d,T\dT2 



(27T) 



-dz- 



sin l j[n - t 2 

r 



(20) 



IV. TUNNELING IN THE PRESENCE OF Senv 

For the purpose of illustration, let us compare S env 
with the action of Refs. 18 ' 30 for a particle moving in a 
dissipative environment: 



Scl = 1/47T / dTidT 2 T} 



(R( n )-R(r 2 )) 2 



(n - t 2 ) 



(21) 



This term results from integrating out the slow de- 
grees of freedom of the environment. The long-time re- 
sponse of these modes reveals itself through a non-local 
in time term in the action. If R(r) describes a parti- 
cle moving with a constant velocity, one can substitute 
[R(t) - R(r')]/(r - t') with the velocity R. Then the 
action Scl reduces to 1/27T / dridT2?7R 2 /2 where the in- 
tegrand is a reminiscent of the Rayleigh's dissipation 
function. The Rayleigh function is used in the Euler- 
Lagrange equations to include dissipation 31 : 



d_ fdC 
dt \dk 



d£ dT_ 
dR + ~dR 



0, 



(22) 



where T = r]R 2 /2 is the Rayleigh function with a fric- 
tion coefficient 77. Since the Rayleigh function enters the 
equation of motion without a time derivative (unlike the 
Lagrangian) , the inclusion of the dissipation into the ac- 
tion costs an additional time integration. Therefore, the 
corresponding term must be non-local in time. 



Under the same approximation of constant velocity, 
the action S env can be written in the coordinate repre- 
sentation as: 



1 

4-7T 



d 2 rdz 



dT X dTi — W 



,(r,z) 



(23) 

The combination \v vor x B„ or (r)] /c is the electric field 
created by a moving vortex. Then, the integrand in 
the action S env is merely the energy dissipation rate 
in the gate caused by the vortex motion, S env = 
\/2n J d 2 rdz J dridT2(jE 2 or (r, z)/2. This expression is in 
full correspondence with the one obtained for a constant 
motion of a particle in the presence of friction rj. 

We return now to Eq. (20) and analyze it for the case 
of tunneling between two minima separated by a distance 
Sr. We start with the integration over the coordinate z 
and momentum q using for the magnetic field B„ or the 
solution given by Eq. (2). We get in result 



S, 



16e 2 A 2 



dr\dT2 



\ In 



(R(n)-R(r 2 )) s 

(n - t 2 ) 2 
8n 2 a 2D d 



(24) 



Tl - T 2 C 



The appearance of the log-factor is very natural if one 
recalls that we integrate the square of the magnetic field 
and it is well known that B vor ~ 1/r at a distance r 
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from the center of the vortex, when r < A. The time de- 
pendence of the logarithm results from the fact that the 
integration over the momenta is limited to the Ohmic 
regime, \oj n \ <C c 2 q/ (2ira 2D ). The time dependence of 
the logarithm is important because in thin superconduct- 
ing films the ratio d/X can be very small. At large time 
differences, which are essential for low temperatures, the 
logarithmic factor in the action becomes hiX/d. 

Following the standard Renormalization Group (RG) 
procedure one gets that the modified tunneling rate Tq 
is: 

r G (T) = T unG (TT tun f E , (25) 

where the exponent Ke is equal to the dimensionless 
dissipation coefficient 

*- < 26 » 

In Eq. (25), T tun is the time of the under-barrier motion 
of the vortex in the process of tunneling in the absence 
of the gate (i.e., the duration of the fast stage of the 
tunneling discussed in the Introduction). The parame- 
ter T t ~Jj acts as the high-energy cutoff because only slow 
excitations that cannot follow adiabatically the tunnel- 
ing particle contribute to the slow stage of the tunnel- 
ing process (it is assumed that T t ~^ >> T). Since only 
the current in the Ohmic regime contributes to S env , 
there is an additional high-energy cutoff ~ c 2 / (4ira 2D d). 
Therefore, in Eq. (25) r t ~^ should be replaced by f t ~„ = 
min{r t ~^, c 2 /(47rcr 2£) (i)}. Little can be said about r tu „ 
as the effective mass of the vortex and the potential 
of the tunneling barrier depend on the specific proper- 
ties of the superconducting film. As to the energy scale 
c 2 / (4Tra 2D d), it is evaluated to be ~ 10 3 K (the resistivity 
of the gate is about 10 fiClcm). 

The temperature enters in Eq. (25) as a low-energy 
cutoff because the excitations with energy smaller than 
T do not contribute to the action. The dependence of 
the tunneling rate on other factors limiting the time of 
response of the environment can be found from the RG 
analysis. 16,32 Apart from the temperature, such factors 
include the typical energy mismatch between the min- 
ima of the vortex potential e and the tunneling rate itself. 
The energy mismatch e in this problem is equivalent to 
a magnetic field in the Kondo problem. Moreover, low- 
frequency components of the Eddy current with frequen- 
cies smaller than Tq cannot develop when the tunneling 
events arc too frequent. This is why the tunneling rate 
determines its own renormalization in a self-consistent 
way. To include the influence of these factors, one should 
substitute T in Eq. (25), by T max = max{T,TQ,e\. 

For Ke > 1 the temperature dependence given by 
Eq. (25) holds for all temperatures down to zero. Thus, 
for e — > the vortex becomes localized at T = 0. In the 
Kondo problem 32 this happens for a ferromagnetic sign of 
the exchange. The localization occurs because the strong 
response of the environment blocks the tunneling. In the 



opposite case, Ke < 1, the tunneling rate remains finite 
when T,e -> 0: 

rG-^tr™^™) 1 ^'. (27) 

This quantity has the meaning of the Kondo temperature 
which determines the physics of switching between two 
states when both T and e are smaller than Tq. Unlike the 
Kondo problem, the dimensionless parameters Ke and 
(^unG T tun) which determine the tunneling rate are com- 
pletely uncorrelated. We are interested in studying the 
case when the tunneling is small even in the absence of 
the gate, i.e., when T uu q <C l/T tun - Then, the tunneling 
between two states is equivalent to the Kondo problem 
in the limit of extreme anisotropy. In this case the renor- 
malization of Ke can be neglected, and the trajectories 
of the RG phase-diagram are straight lines. 

So far, in the discussion of the renormalization of the 
tunneling rate we ignored the time dependence of the log- 
factor in Eq. (24). Generally speaking, this is not valid as 
there can be a big window of energies, f 4 ~„ > \u> n \ > t7 , 
where t^ 1 = c 2 /Ana 20 X. Then, the RG analysis has to 
be revised. The modified phase diagram is plotted in 
Fig. 3; see Appendix for details. Notice, that the con- 
dition for localization (the inability to tunnel) becomes 
harder. It is also worth noting that the Delocalization- 
Localization phase transition occurs at a finite energy 
scale. 

r t 

A unG v turu 







derealization / 




/ 




localization 











FIG. 3: The Delocalization-Localization phase diagram for the 
tunneling of vortices in the plane of T un GT t un and the dissipa- 
tion coefficient Ke- 



Let us discuss the result obtained for the exponent Ke- 
Naturally, the action S env is proportional to the square 
of the magnetic field. After integration over the coordi- 
nates one obtains the factor (a/ A) 2 . From the structure 
of the action (24) it follows that Ke is proportional to 
the squared distance of the change of vortex position, 
which is typically of order £. Altogether, the exponent 
Ke acquires the factor (£/A) 2 <C 1. This smallness is op- 
posed by the sheet conductance of the gate, a 2D , which 
in the case of a thick well conducting gate gives a large 
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factor 2irh<jd/e 2 . As usually, the effect of dissipation is 
stronger for cleaner systems. 



V. ORTHOGONALITY CATASTROPHE 

The magnetic flux of a vortex piercing through the gate 
scatters the electrons in a way similar to Aharonov-Bohm 
(A-B) scattering. 21 ' 22 The tunneling of a vortex enforces 
the electrons inside the gate to adjust to its new posi- 
tion. The response of the electrons to a sudden change 
of the vortex position leads to the OC that manifests 
itself in the vanishing overlap of the two wave 

functions describing the macroscopic electron system be- 
fore and after the change of the scattering potential. 19 
The tunneling rate renormalized by the overlap integral 
is 



r 



G 



unG 



(28) 



The overlap (^/l^i) can be expressed in terms of the 
operators Si and Sf describing the scattering of the elec- 
trons by the magnetic field of the vortex in its initial and 
final positions 33 : 



\(*f\*i)\=N- K ° ; 

K oc = -^Tr{\n 2 (S f S^)) 



(29) 



Here N is the number of electrons in the gate and, 
hence, the overlap factor vanishes unless there is a mech- 
anism that limits the effectiveness of the OC. It is clear 
from the comments to Eq. (25) that at finite tempera- 
tures 34 ' 35 the parameter 1/N should be substituted by 
(max{T, Ta^tun)- Obviously, the renormalized tunnel- 
ing rate Tq is given by Eqs. (25) and (27) with Kb re- 
placed by Kqci an d the localization of vortices can be 
achieved when Koc > 1- 

In the following part of this section we show that for a 
superconducting film magnetically coupled to a metallic 
gate (see Fig. 2) the exponent Koc is 



K. 



oc 



gate \ 
F I 



(a5r) 2 
64A 



(30) 



where k 9 F ate is the Fermi momentum of the electrons in 
the gate and the prefactor <; is evaluated numerically as 
« 0.4. 

The cylindrical symmetry of the vortex allows us to an- 
alyze the scattering of electrons using the basis of cylin- 
drical waves, \£,q,k z ); here £ is the angular momentum 
along the z axis, while q and k z are the magnitudes of 
the in-plane and z components of the momentum. In 
this basis the elements of the matrix SfS~ l can be easily 
calculated in terms of the phase shifts Si as 



f {e\s f srY)f 



(31) 



-2iS„ 



' : J n {q6r)J n ^ e , +e (q6r) 



where J v (z) is the Bessel function. 

To proceed further, we need to find the specific phase 
shifts for scattering by a vortex. An analogy to classical 
scattering, where the angular momentum is related to the 
impact parameter b = \£\/q, helps elucidate the behavior 
of the phase shift as a function of £. For b ^> A, the scat- 
tering by the vortex is similar to the A-B scattering by 
a flux a$o. In the A-B scattering 21,22 electrons acquire 
the phase 8f~ B = ^{\£\ — \£ — a\). The uniqueness of 
this scattering is in its infinite range: 8% does not vanish 
when \£\ — > oo. For scattering by the vortex, the jump 
in the A-B phase shifts is smeared out, but the infinite 
range character of this scattering is preserved. Hence, 5e 
varies monotonically as a function of £ between the two 
limits: 



8e ► a— son £. 

t»q\ 2 



(32) 



Naturally, for qX 3> 1 the phase shift depends on b and A 
only through the dimensionless combination b/X = IjqX 
such that 5 t = ^g(£/qX); see Fig. 3 in Ref. 20 for 
illustration. 

We now notice that the sum determining the elements 
of f{£\SfS~ 1 \£')f is accumulated at — qSr < n < qSr. 
This is because the Bessel Functions J u {z) decay expo- 
nentially with their order when v > z. Therefore, since 
Sr/X ~ £/A <C 1, the phase shifts difference in Eq. (31) 
can be approximated as: 



Jn+e 



r/ r/ a7r i ( £ 



< 1. (33) 



The final step of the calculation is to expand in Sr/X 
the logarithm in Eq. (29), and take the trace over £ and 
the momentum on the Fermi surface. The outcome of 
the calculation is given in Eq. (30). The gate thickness d 
appears here as a result of taking the trace. The specifics 
of the vortex solution enter only through g(x), with the 
integral yielding q = J dx(dg/dx) w 0.4. 

Using the known expressions for A and £ in disordered 
thin films, the exponent can be rewritten as: 



Kqc ~ C 



48tt 



(k 9 F ate d)(k 9 F ate a)(k s F c r) 



(34) 



The index sc refers to the electrons in the superconduct- 
ing film: l sc is their mean free path (in the normal state) 
and Dp is the Fermi velocity. Interestingly, T c drops out 
from Koc so that it depends only on the geometrical fac- 
tors and the non-superconducting properties of electrons. 
We see that the value of the exponent Koc is determined 
by a small factor ~ 10~ 7 opposed by a product of a few 
large factors. Unlike Ke, the condition for vortex local- 
ization can be easily fulfilled by Koc for a not too thin 
gate and not too disordered superconducting film. 
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VI. COMPARISON BETWEEN THE TWO 
CALCULATIONS AND DISCUSSION 

The expressions in Eqs. (26) and (30) for the exponent 
describing the rcnormalization of the tunneling rate have 
been obtained assuming that only one vortex participates 
in each tunneling event. In general, vortices can tunnel 
as a bundle, or as topological defects. In a vortex lattice 
or in a glass state such defects can be dislocations pairs, 
interstitials or vacancies. Still, the calculation remains 
valid as long as 5r <C A. This is because the magnetic 
field of the tunneling vortices extends over a large dis- 
tance, so that their exact configurations before and after 
the tunneling are not important. The only relevant quan- 
tity is the product adr. For a single vortex, a = 1/2. 
When more than one vortex tunnel together a should be 
multiplied by the number of vortices. 

The two expressions for the exponent, Ke and Koc, 
share the same dependence on adr and d. Therefore, the 
ratio between them is " universal" : 

K E crln(X/d) I 

qate /or\ 

Ko~c~ ~ e 2 (k g F at y\ ~ ~T' 1 ' 

where l ga te is the mean free path of the electrons in the 
gate. Since in a thin superconducting film the penetra- 
tion depth A is very large, under the conditions of the 
experiment 9 , the OC is dominant. 

In order to understand the difference between the two 
expressions for the exponent, we have to explain the de- 
pendence of Koc on A. Although we invoke the expan- 
sion in terms of the small parameter Sr/X <C 1, we get 
Koc (Sr) 2 /X. This is typical for the OC when an ex- 
tended scattering potential is considered, because in this 
case a large number of scattering channels (harmonics) is 
involved. Therefore, the total effect of OC is paramctri- 
cally bigger than the one from a single channel. 36 For the 
problem discussed here, the relatively weak dependence 
of Koc on A can be understood from the following ar- 
guments. It has been shown that the OC is determined 
by X^(<^ — <^+i) 2 ~ J2e(fi'e) 2 - Since the phase shifts 
approach asymptotically the limit ia^, the sum 

Therefore, the result obtained for the exponent Koc cor- 
responds to the differences (Se — Se+i) that are distributed 
almost equally between L ~ qX channels: 

E^-<W~£(f) 2 ~y. 07) 

Indeed we see that the first power of A -1 is natural for 
the exponent Koc because a large number of scattering 
channels is perturbed by the change in the potential when 
such an extended object as a vortex tunnels. 

The peculiarity of the discussed problem is that a vor- 
tex in a thin film is a very extended object. In gen- 
eral, the tunneling of an extended object excites many 



channels of the environment. Unlike the standard OC 
caused by a scattering potential (but not a vector poten- 
tial) where the exponent is determined by sum of squares 
of the phase shifts, here the sum of the squared deriva- 
tives of the phase shifts determines the exponent of the 
OC. Still there is some similarity between the two prob- 
lems. In the standard OC the sum of the phase shifts is 
finite because of the Friedel sum rule (see e.g. the tun- 
neling via the localized level considered in Rcf . 36) , while 
in the discussed problem the asymptotic limits of the AB 
phase shifts make the sum of the derivatives of the phase 
shifts to be finite, see Eq. (36). As a result of this "AB 
sum rule" the OC exponent given by Eq. (37) is similar 
to that given by Eq. (10) in Ref. 36. 

One may conclude from Eqs. (36) and (37) that the 
randomization of the phase shifts due to the disorder 
can only increase the value of the exponent Koc- (In 
the general case, £ should be substituted by the index 
of the states diagonalizing the scattering matrix.) The 
scattering by impurities leads to the randomization of 
the phase differences, while the asymptotic limits of the 
phase shifts remain the same, ±ck^. Therefore, the value 
of the exponent Koc, which is determined by the squares 
of the phase differences, should increase in the presence 
of disorder. This conclusion is in accordance with the 
existing theoretical results on the enhancement of the 
OC-exponent by not too strong disorder. 37 ' 38 

The sum over the channels (harmonics) enters in a nat- 
ural way into the exponent Koc, while in the calculation 
of the "tunneling with dissipation" the sum is absent. 
The scheme of calculation of the "tunneling with dissi- 
pation" for the Eddy currents corresponds to the OC 
expression calculated in the perturbation theory up to 
the second order with respect to the change of the vector 
potential. It is merely finding the elements of the scatter- 
ing matrix in the Born approximation. However, the AB 
effect is a non-perturbative phenomenon. When the cur- 
rent is calculated in the linear response the asymptotic 
limits of the phase shifts are not captured. One may also 
see that in the linear response scheme only one channel 
remains in the sum over the channels that is present in 
the OC expression. Therefore, the additional factor of A 
cannot be reproduced within the "tunneling with dissipa- 
tion" scheme. This is the reason why Koc is significantly 
larger than the obtained Ke- 

The idea to use a double layer system to study the dy- 
namics of vortices is well known. 39-41 In addition to the 
magnetic coupling between the film and the gate, one 
may consider a capacitive coupling between them. In the 
case of the Joscphson junction arrays (or granular super- 
conductors) the capacitive coupling reduces the fluctua- 
tions of the phase of the superconducting order parame- 
ter. 41 ' 42 As a result, the system may undergo a transition 
from an insulating to a superconducting state. However, 
for a homogenous film with a relatively small resistance 
~ 1.5fcf2/D used in Ref. 9 the phase fluctuations are not 
so effective. 43 This is confirmed by the observed insensi- 
tivity of the critical magnetic field H c to the presence of 
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the gate. Furthermore, in homogeneous superconductors 
the motion of vortices is not accompanied by the redistri- 
bution of the charge density. Therefore, there are good 
reasons to ignore here the capacitive coupling between 
the film and the gate. 



VII. SUMMARY 

In this work we have studied how a metallic gate placed 
above a superconducting film affects the tunneling rate 
of the vortices. The gate and the film are coupled by 
the magnetic field of the vortices that pierces through 
the gate. We analyze the renormalization of the tun- 
neling rate by the gate. We consider two approaches to 
describe the response of the electrons inside the gate on 
the tunneling event: (i) the Eddy current in the gate gen- 
erated by the motion of vortices, and (ii) the OC caused 
by the change in the vortex position. The OC is due 
to the Aharonov-Bohm scattering of electrons inside the 
gate. The exponent determining the renormalized tun- 
neling rate T G (T) = T unG (Tr tM „) A ' is given by Eq. (26) 
for the effect of the Eddy current, and by Eq. (30) for the 
OC. We find that for the experimental setup of Ref. 9 the 
effect of the OC provides an exponent sufficient for a sub- 
stantial suppression of the tunneling rate of the vortices, 
K oc ~ 1. 

The peculiarity of the discussed problem results from 
combination of two elements: the extended size of the 
tunneling vortex and the unique features of the AB scat- 
tering. In general, the tunneling of an extended object 
creates excitations in many channels of the environment. 
The sum over the channels (harmonics) enters in a nat- 
ural way into the exponent describing the effect of the 
OC, while in the calculation of the "tunneling with dissi- 
pation" the sum is absent. This is the reason why Koc 
is significantly larger than the obtained Ke- The scheme 
of calculation of the " tunneling with dissipation" for the 
Eddy currents corresponds to the OC expression calcu- 
lated in the perturbation theory up to the second order 
with respect to the change of the vector potential. Be- 
cause of the non-perturbativc character of the AB effect, 
the phase shifts cannot be found within the Born approx- 
imation. Therefore the action S env , being formulated in 
terms of the macroscopic charge and current densities in 
the regime of the linear response, is unable to describe 
the response of all fluctuation modes that can be excited 
by a vortex as a result of tunneling. 

We address our analysis to a recent experiment in 
which the resistance of a superconducting film has been 
measured in a magnetic field both with and without a 
gate. We interpret the experiment by assuming that the 
origin of the resistance at low temperatures is tunneling 
vortices. From our point of view, the difference in the 
resistance of the gated and ungated film indicates that 
the gate reduces the tunneling rate of the vortices mak- 
ing them localized. Indeed, we show here that adding a 
gate may effectively suppress the vortex tunneling. The 



gated system discussed here can be used as an effective 
experimental tool for investigating the vortex motion at 
low temperatures. The gated system provides a unique 
opportunity to study the vortex tunneling in thin super- 
conducting films by such simple means as varying the 
characteristics of the gate, in particular the gate thick- 
ness and/or the sheet conductance of the gate. This may 
help to identify the different mechanisms that contribute 
to the suppression of the vortex tunneling rate. 
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APPENDIX A: 
LOCALIZATION-DELOCALIZATION PHASE 
DIAGRAM FOR EQ. (24) 

The renormalized tunneling rate as given by Eq. (25) 
corresponds to the RG equation 

d\n{T G T)/d\n{T/f tun ) = (1 - K E ). (Al) 

This equation is valid as long as the logarithmic factor 
in Eq. (24) is independent of the time difference \t\ — t?\. 
The time dependence in the logarithmic factor is essen- 
tial for the energy interval r t ~„ > 1/t > t^ 1 , where 
t^ 1 = c? I AiradX. Then, the RG equation has to be 
modified by replacing the dissipation coefficient Ke in 
the above equation by the energy dependent parameter 
K E = K E [l-\n(Tx/ft un )/\n(X/d) + g/\n(X/d)], with the 
logarithmic variable g = ln(r/rt un ). For the discussed 
energy interval the "RG equation" becomes 

dMT G r)/dg=(l-k E {g)) (A2) 

One should not be confused with the appearance of g 
in the RHS of the "RG equation". Here we just inte- 
grate S env in the exponent determining the renormalized 
tunneling rate T G oc exp(— S env ). In a sense we calcu- 
late something like the Dcbye- Waller factor created by 
the dissipative environment 44 . [The g-dependence of the 
parameter K E can be reformulated as an RG equation, 
additional to Eq. (A2). Namely, 

dk E /dg = { c T l ° < f f h ; (n/ f t[m); (A2a) 

where const = K E /ln(X/d). Because of the different 
dependencies of Ke on g, the RG process may develop 
in two steps. We discuss the details below.] 

Localization occurs if in the course of the RG pro- 
cess the tunneling rate T G decreases faster than the run- 
ning scale r -1 . Hence, the line separating the localized 
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and dclocalizcd states is determined by the condition 
Ke(J"g) = 1 as long as T G > t^ 1 . Since the expres- 
sion for Ke assumes that g < \n(T\/ft un ), the modified 
Ke is smaller than Ke- Therefore, the condition for lo- 
calization (the inability to tunnel) becomes harder. It 
is worth noting that in the discussed problem the phase 
transition occurs at a finite energy scale 

ln(l/r G ft tt „) = l I^± _ [ln(A/d) -ln(r A /f t „„)]. (A3) 
Ke 

To avoid unnecessary complications, from now on we 
limit ourselves to the case when the energy cutoff 1 jftun 
is determined by c 2 / Anad 2 . Then ln(T\/f tun ) = \n(X/d) 
andEq. (A3) reduces to ln(l/r G f tlm ) = In (T\/f tun )/K E - 
Note that the line of the phase transition exists for all 
Ke > 1 and in the delocalized phase r G > t^ 1 . Clearly 
for K E < 1 the tunneling rate always remains finite. 

In order to find the line of the Delocalization- 
Localization transition in the plane of the dimcnsionlcss 
parameters (TunGTtun) and Ke, one has to integrate back 
Eq. (A2) starting at ln(r G r) = 0. (The value of T G cor- 
responding to the line of transition should be found from 
the condition K e (Tg) = 1)- This procedure yields for 
the boundary between the two phases 

ln(l/r ) = \n(T X /f tun )/2K E . (A4) 

One may see that along the boundary (r G r tu „) = 

(r ^ unG^tun) ■ 

In the delocalized phase, the renormalized tunneling 
rate can be found by integrating the RG equation starting 
from the bare tunneling rate down to the energy when 
ln(F G r) = 0. Since for K E > 1 the condition ln(r G r) = 
is satisfied before the running scale r _1 reaches tF , 
the RG process involves only Eq. (A2). For Ke < 1 the 



situation is more delicate. In the first step of the RG 
process Eq. (A2) is integrated. If r _1 reaches tF 1 before 
ln(r G r) = (i.e., T G < rF 1 ), the process should be 
continued. (In terms of the bare parameters, this occurs 
when r imG < T^ l {T\/f tU n) KE ^ 2 )- In the second step one 
has to integrate Eq. (Al) using t^ 1 as an upper cutoff 
instead of f t ~ n and r G (r^" 1 ) = T\ as an initial value; here 
T\ = T un G(ftun/i~\) KE t 2 is the result of the integration 
in the previous step for T G . Finally, after the two-steps 
rcnormalizion the tunneling rate is 

r G ~ tf 1 (r A T A ) 1/(1 -^ } , (A5) 

which is equivalent to 

r G f tun ~ (r unG T tun ) 1 ^ 1 - K ^ {fx/f tun ) KE/2{1 ' KE) . 

(A6) 

The Delocalization-Localization phase diagram for the 
tunneling of vortices in the plane of T un QT tun and the dis- 
sipation coefficient Ke at zero temperature is presented 
in Fig. 3. 

Finally, we touch upon the role of the relativistic ef- 
fects. In the derivation of the effective action (24) the 
integration over the momenta has been limited to the 
Ohmic regime, q > 2wa 2D \uj n \/c 2 . We are interested in 
a gate with high conductivity, such that 2ir a 2D /c 3> 1. 
In this case there are relativistic effects that have been 
left out from the action. As it has been mentioned in 
Sec. Ill, at high frequencies |w„|/c > q the dissipation 
occurs through the Cherenkov radiation. The produced 
dissipation is equivalent to an effective two-dimensional 
conductivity that is approximately c/2tt. Therefore, this 
source of the dissipation is negligible compared to the 
Ohmic dissipation and cannot change our conclusions. 
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